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Abstract. 

In the coming decade, the LIGO/VIRGO network of ground-based kilometer- 
scale laser interferometer gravitational wave detectors will open up a new 
Astronomical window on the Universe: gravitational waves in the frequency 
band 10 to 10 4 Hz. In addition, if the proposed, 5 million kilometer long, 
space based interferometer LISA flies, another window will be opened in 
the frequency band 10 -4 to 1 Hz. We review the various possible sources 
that might be detected in these frequency bands, and the information that 
might be obtainable from observed sources. Several key possible sources 
are inspirals and coalescences of neutron-star neutron-star and/or neutron- 
star black-hole binaries; inspirals, mergers, and ringdowns of black-hole 
black-hole binaries (both solar mass and supermassive); stellar core col- 
lapse; rapidly rotating neutron stars; the formation of supermassive black 
holes; and inspirals of compact objects into supermassive black holes. 



1. Introduction 

This review of gravitational wave sources is divided into three sections. 
First, we review the detector sensitivities that have been achieved to date 
and discuss projected sensitivities for detectors now under construction. 
Second, we summarize the current observational upper limits on gravita- 
tional waves in various frequency bands. 

The main body of this review will consist of a survey of various antic- 
ipated sources of waves. Each anticipated source can be roughly charac- 
terized by a characteristic frequency / and a characteristic value of strain 
amplitude h. However, it is important to also note that sources vary widely 
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with respect to how uncertain is the rate of their occurrence in the Uni- 
verse. The enterprise of anticipating potential gravitational wave sources is 
extremely uncertain. For most sources that we can conceive of, either the 
wave strengths are uncertain by several orders of magnitude, or the event 
rate is uncertain by several orders of magnitude, or the very existence of the 
source itself is very uncertain. While there are some important exceptions 
such as coalescing compact binaries, these are the exception rather than the 
rule. The upside of this great uncertainty is the potential for gravitational 
wave astronomy to bring us new and interesting information. 

For more details on the topics discussed here, the reader is encouraged 
to consult the detailed recent review articles by Thorne [1, 2, 3, 4], and also 
the review article by Allen [5] on stochastic gravitational waves. 

2. Detector Sensitivities: Achieved and Projected 

As is well known, there are two technologies which are being pursued for 
detection of gravitational waves. The first is to monitor the modes of vi- 
bration of a solid test mass which is cooled to very low temperatures. Such 
detectors are called resonant mass detectors, or more colloquially "bars". 
The second technology is to monitor the relative displacement of suspended 
test masses using laser interferometry. 

There are presently operating several cylindrical resonant mass anten- 
nae: the ALLEGRO detector at Louisiana State University [6], the NIOBE 
detector at the University of Perth, Australia [7], the EXPLORER de- 
tector at the University of Rome [8]. All of these detectors are operating 
at temperatures of a few degrees Kelvin and have sensitivities of order 
h ~ (several) x 10~ 19 . There are also one or two so-called third gener- 
ation detectors which are operating at temperatures of a few tens of mi- 
liKelvin: the AURIGA detector at the University of Legarno, Italy [9], and 
the NAUTILUS detector at the University of Rome [10]. These detectors 
have sensitivities of order h ~ (several) x 10~ 20 [10]. 

In addition, there are plans to construct a new generation of resonant 
mass detectors which will be spherical, or nearly spherical (TIGAs or Trun- 
cated Icosahedral Gravitational Arrays [12]), or hollow spheres. Detectors 
in the planning stage include the GRAIL project in the Netherlands [13], 
the OMNI-1 project in Brazil, and the TIGA project in the USA [12]. 
These detectors have design sensitivities of order 10~ 21 ; see Fig. 1 below. 
Attaining this sensitivity level will require reducing the so-called noise tem- 
perature by several orders of magnitude, and so will not be an easy task. 
More details on resonant mass detectors can be found in the contribution 
by Massimo Cerdonio to this proceedings. 

Interferometer detectors have several advantages over resonant mass de- 
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tectors, the primary one being that they are intrinsically broadband while 
resonant detectors are narrowband. The technology for such detectors has 
been continuously under development for several decades. Today an in- 
ternational network of kilometer-scale, interferometer detectors is under 
construction, consisting of: 

• The American LIGO project [11], which initially will consist of 
three separate interferometers, an interferometer with 4 km arms 
in Livingston, Louisiana, and two interferometers with armlengths 
of 2 km and 4 km in Hanford, Washington. The detectors will come 
online around 2001. 

• The French- Italian VIRGO project [14], which is constructing an 
interferometer with 3 km arms near Pisa, Italy. This detector is 
also planned to come online around 2001. 

• The British-German GEO600 project [15], which is constructing 
a 600 meter long interferometer near Hanover, Germany, coming 
online around the year 2000. 

• The Japanese TAMA project [16], which is constructing a 300 
meter interferometer near Tokyo. 

There are also several prototype interferometers, for instance, the LIGO 
40m prototype and the Garching and Glasgow interferometers associated 
with the GEO 600 project. The kilometer-scale interferometers will act as a 
coordinated network: for each burst of waves, the data from all the detectors 
will be analyzed to obtain the two independent waveforms h + (t) and h x (t), 
and also the direction to the source. More details on this network can be 
found in the contribution by Noma Robertson to these proceedings. 

Both resonant mass detectors and earth-based interferometer detec- 
tors are limited to the frequency band 10 Hz <, f <> 10 4 Hz. At lower 
frequencies, the background of time varying near-zone Newtonian gravi- 
tational fields produced by objects on the Earths surface (winds in the 
atmosphere, seismic waves in the Earth's crust etc.) becomes very large. In 
order to attempt to measure gravitational waves in the band of frequencies 
10~ 4 Hz <, f <> 1 Hz, it is necessary to go into space. Laser interferometry 
can be used to monitor the distances between spacecraft and thus to mea- 
sure gravitational waves. One proposed mission of this kind is the Laser 
Interferometer Space Antenna (LISA), which the European Space Agency 
tentatively plans to fly around 2015 [17]. LISA consists of six spacecraft in 
a triangular configuration in a solar orbit. 

Figure 1 shows the noise levels that have been achieved to date in several 
representative detectors: bar detectors in 1976, in 1990 (Louisiana/Rome), 
in 1997, and the Caltech 40 meter prototype interferometer in 1990 and 
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Figure 1. This plot shows the noise levels that have been achieved to date in several 
gravitational wave detector systems: the Caltech 40m interferometer [18], and the best 
bar detectors in 1976, [1], in the early 1990's [2], and in 1997 [10]. It also shows that 
projected noise levels for future detector systems: initial and advanced LIGO interfer- 
ometers [11], the spacebased LISA interferometer [2], and future generations of resonant 
mass detectors in the form of truncated icosahedral gravitational arrays (TIGAs) [20]. 
For the broadband interferometer detectors the quantity plotted is /i m8 (/) = -J fSh(f)- 
For the narrowband resonant mass detectors (bars and TIGAs), the quantity plotted is 

ftrms.bar = a/irms (/o) \J /o/A/, where A/ is the bandwidth of the noise spectrum, fo is 
the central frequency, and a is a constant of order unity chosen so that narrowband and 
broadband detectors will have the same value of h on this plot if they are equally efficient 
at detecting broadband, burst waves [21]. 



1994. Also shown in Fig. 1 are the projected noise levels for future detec- 
tors: the first generation "initial" LIGO interferometers, later "advanced" 
LIGO interferometers, the LISA interferometer, and proposed TIGA reso- 
nant mass detectors. There will also be a generation of interferometers in 
LIGO intermediate between the initial interferometers and the advanced 
interferometers; these "enhanced interferometers" will be obtained by in- 
crementally improving the various elements of the initial interferometers 
over a time scale of a few years. In Fig. 2 we compare the noise spectra for 
all the major interferometer projects. 

For the broadband, interferometer detectors, the quantity plotted in 
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Figs. 1 and 2 is the the rms dimensionless strain per unit logarithmic fre- 
quency, forms (/)• This is defined so that the total noise-induced fluctuation 
in the output h(t) of any detector can be written as 

/oo 
d(ln/)forms(/) 2 . (1) 
-oo 

For the resonant mass detectors, on the other hand, the quantity l/h rms (f) 
sharply peaked, and the signal-to-noise ratio for any broadband burst of 
waves is proportional to f dfl/h Yms (f) 2 = Af /fo rm s(/o) 2 > where Af is the 
bandwidth and fo is the central frequency. In Fig. 1 we plot for resonant 
mass detectors the quantity 




forms, bar — \\ » r ^rms(/o)i (2) 



where a is a constant of order unity chosen so that narrowband and broad- 
band detectors will have the same value of fo if they are equally efficient 
at detecting broadband, burst waves [21]. [Note that, when considering the 
detection of periodic signals, one should instead compare directly the values 
of forms(/) °f bars and interferometers.] 

For a given broadband burst of waves, broadband detectors can measure 
the gravitational waveform h(t), whereas narrowband detectors can essen- 
tially measure only the Fourier transform of the waveform at the resonant 
frequency. Although we have treated the planned TIG A/spherical detec- 
tors as narrowband instruments above, in fact they will have Af / fo ~ 0.1 
[12], in contrast to bar detectors to date which have A///o ~ 10~ 3 . A 
"xylophone" of TIGAs could thus act as a broadband instrument [12]. 

Consider now how strong waves have to be in order to be detected. There 
are three different types of gravitational waves: burst waves, periodic waves 
and stochastic waves. For each type of wave, let h amp (t) denote the physical 
size of the metric perturbation. For burst waves, a signal with characteristic 
frequency f c will be detectable whenever fo c fr ^ h Tms (f c ), where the effective 
strain amplitude fo e ff is 



^ ff "^ m V21n[r/At]- (3) 

Here A/" cy cies is the number of cycles of the waveform in the bandwidth of 
the detector, T is the total duration of the data set in which one searches 
for a signal, and At is the effective sampling time [23]. For periodic waves, 
a signal with frequency f c will be detectable when h e s <i h ITas (f c ), where 



VTAf 
a/2111 [TAf] 



foeff = foamp— 7= \m \ ¥\ (^) 
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Noise Power Spectra 
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Figure 2. The projected noise spectra for the current major interferometer projects. 
The LIGO spectra were taken from Fig. 10 of Ref. [11] and the spectra for the other 
interferometers were obtained from Ref. [22]. 

and A/ is the range of frequencies one searches through. Finally, stochastic 
waves will be approximately detectable when h e s h rras , where 



and Af is the bandwidth over which the detector is sensitive. Plotting 
h e g versus h rms for several different detectors and various different sources 
of waves illustrates the waves' detectability. The definitions (3) - (5) are 
approximate but correct to within factors of order unity; more precise char- 
acterizations of detectability for burst, periodic and stochastic waves can 
be found in Refs. [1], [70] and [5], respectively. 

3. Indirect detections and observational upper limits 

The most celebrated indirect detection of gravitational waves is of course 
the monitoring of the orbital decay of the binary pulsar PSR 1913+16, for 
which the decay rate agrees with the prediction of general relativity to bet- 
ter than 1% [24]. Aside from this, several different techniques have been 
used to place upper limits on the energy density in gravitational waves in 
different frequency bands. One technique is to monitoring the radio sig- 
nals from interplanetary spacecraft; such signals would be Doppler shifted 
in a characteristic way by a burst of gravitational waves. The limits ob- 
tained from this technique [1] are illustrated in Fig. 3. Another technique 



h c S — fr-amp \/T/S.f 



(5) 
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Figure 3. Current observational upper limits on gravitational wave energy density in 
various frequency bands, together with the projected sensitivities of advanced LIGO 
interferometers and of LISA. 



is to monitor precisely the arrival times of radio pulses from pulsars (ac- 
curacy ~ 1 fis) over several years; the resulting limits [25] are shown in 
Fig. 3. From primordial nucleosynthesis we obtain the constraint on pri- 
mordial gravitational waves that J din f £l g (f) ^ 10 -5 [5]; this limit does 
not apply to waves that were produced after neucleosynthesis. Finally, it is 
likely that some portion of the anisotropy of the cosmic microwave back- 
ground radiation is due to primordial gravitational waves (the remaining 
portion being due to density or scalar perturbations). The current obser- 
vations of the anisotropy spectrum cannot distinguish between scalar and 
gravitational wave contributions. Thus, these observations provide upper 
bounds on £l g (f) at low frequencies as illustrated in Fig. 3 [5]. They may 
also constitute an indirect detection of gravitational waves, although cur- 
rent theoretical prejudice suggests that the gravitational wave contribution 
to the anisotropy is smaller than the noise level in the observations to date. 

The limits in Fig. 3 are shown in terms of the quantity which 
is the energy density in gravitational waves per unit logarithmic frequency 
divided by the critical energy density /^critical necessary to close the Universe: 

1 dW 

n g (f) = (/) 

/Ocritical d\nfd 6 X 

( f \ 2 ( fram P \ 2 ( 100km s-iMpc- 
UOOHzy U-3xl0- 20 / ^ H 
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where Hq is the Hubble constant. 

4. Overview of Sources 

As discussed above, the the high frequency band of 10 Hz <> f <^ 10 4 Hz is 
the domain of earth-based instruments. Key sources in this band are stel- 
lar core collapses, coalescences of compact binaries (neutron star - neutron 
star, neutron star - black hole, and black hole - black hole binaries), vi- 
brating/precessing neutron stars, and black hole births. The low frequency 
band of 10~ 4 Hz £ / £ 1CT 1 Hz is the domain of LISA. Here the sources 
are binary stars, formation and coalescences of supermassive black holes 
(SMBHs), and capture of compact objects by SMBHs. Sources of stochas- 
tic waves include parametric amplification of vacuum fluctuations during 
inflation; defects such as cosmic strings, and phase transitions in the early 
Universe. These sources can produce waves over a huge frequency band 
from / ~ 10~ 18 Hz to / ~ MHz. In addition, the superposition of very 
many indistinguishable coalescing compact binary signals or stellar core 
collapse signals may produce a detectable stochastic background in the high 
frequency band [26]. Finally, more speculative sources include theoretical 
constructs such as naked singularities, boson stars, soliton stars, etc. 

5. Stellar Core Collapse 

Type II supernovae occur roughly once every 40 years in our Galaxy. The 
rate of stellar core collapses could be somewhat larger than this, since some 
stellar core collapses could be electromagnetically quiet and/or associated 
with the accretion induced collapse of white dwarfs. While the electro- 
magnetic signal is dominated by the ejected mantle, the gravitational wave 
signal is dominated by the dynamics of the collapsing core. Numerical mod- 
eling of the dynamics of core collapse and bounce is complex: despite much 
work in this area there is not yet a firm consensus. The amount and charac- 
teristics of the gravitational wave emission depends on hydrodynamic pro- 
cesses that have not yet been well modeled, and also on the initial rotation 
rate of the degenerate stellar core before collapse, which is poorly known. 
Thus, there are large uncertainties associated with this type of source. In 
fact, our knowledge is so poor that the event detection rate could be as 
large as many per year with initial LIGO interferometers, or as poor as less 
than one per year with advanced LIGO interferometers. 

There is a critical value of the initial angular momentum J of the degen- 
erate stellar core. If J = 0, then the core will bounce once nuclear densities 
are reached at r ~ 10 km, where r is the radius. Otherwise the eccentricity 
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e increases during the collapse; to a good approximation 

T 2 J 2 1 



(7) 



where T is the rotational kinetic energy of the core, W is its gravitational 
potential energy, J is its conserved angular momentum, and M is its mass. 
The bounce occurs at eccentricities of order unity when 



which corresponds to initial rotation periods of a few seconds. For J 3> 
J cr it, centrifugal flattening will hangup the collapse. On the other hand, 
if J <C J C rit) there is no centrifugal hangup and it is thought that core 
remains axisymmetric during the collapse. Numerical simulations predict 
that in this case the total energy radiated from the collapse and bounce is 
A-Erad ~ W~ 7 M Q c 2 , and the effective strain amplitude is [27] 



Such a core collapse would be visible only in the local group of galaxies by 



Consider now the case in which J 3> J cr it- Here centrifugal flattening 
and hangup occur before nuclear densities are reached; the protoneutron 
star could form a disc with radius ~ 100 km or larger. If T/W ^ 0.3, 
theory suggests that a the fluid configuration is unstable to becoming non- 
axisymmetric via a dynamical triaxial instability. This would greatly en- 
hance the gravitational wave emission; a purely axisymmetric body can 
emit waves only via its collapsing motion, whereas a non-axisymmetric 
body can emit waves via its rotational motion. Such a triaxial instability 
has been seen in numerical simulations [28], in which the hydrodynamical 
processes are treated in a simple and approximate way. [Note however that 
one recent numerical simulation [29] finds that the non-axisymmetry does 
not appreciably enhance the gravitational wave emission contrary to expec- 
tations.] If such an instability occurs, as much as ~ W~ 3 M Q c 2 of energy 
could be radiated into gravitational waves, and such collapses could be seen 
out to the VIRGO cluster by initial LIGO interferometers and to several 
hundred Mpc by advanced LIGO interferometers. However, the fraction of 
core collapses (if any) which become very non-axisymmetric is unknown. 
Nevertheless, even if only one core collapse in 10 4 undergoes a tri-axial in- 
stability, such collapses might be the most common type of core collapse 
seen by the detectors. Several specific scenarios for non-axisymmetry are 
reviewed by Thorne [4]. 
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There is little observational evidence concerning the initial rotation rates 
of degenerate stellar cores. Two points are worth noting, however. First, the 
increasing evidence that neutron stars undergo substantial kicks in super- 
novae [30] (achieving velocities ^> 1500 kms -1 in some cases) suggests that 
core collapses are very non-symmetric. Second, observations suggest that 
newborn neutron stars are typically not born spinning near the breakup 
angular velocity that one would expect if the core collapse has an excess of 
angular momentum [31]. The traditional view has been that this indicates 
that most pre-collapse stellar cores are slowly rotating. However, it has re- 
cently been predicted that newborn neutron stars can spin down from near 
maximal rotation down to ~ 20 Hz within a year or so due to r-mode insta- 
bilities (see Sec. 7 below). Thus, there need not be a contradiction between 
the observational evidence that rapidly spinning (periods ~ milliseconds) 
neutron stars have typically accreted most of their angular momentum, 
and the supposition that most pre-collapse degenerate stellar cores could 
be rapidly rotating. 

Detecting stellar core collapses would have several payoffs for astron- 
omy [1]. If the core collapse were close enough to produce a detectable flux 
of neutrinos, one could constrain neutrino masses by comparing the arrival 
times of gravitational waves and neutrinos. For example, a coincident de- 
tection of gravitational waves and neutrinos from the supernova SN1987a 
would have yielded the limit vn v <> 4eV. One could also obtain informa- 
tion from the gravitational waveforms about the dynamics of the core col- 
lapse and bounce, tipoff optical astronomers before optical supernova first 
brighten (so they can measure the early portion of the light curve), and 
obtain observational evidence about the rate of formation of neutron stars. 

6. Coalescences of Compact Binaries 

6.1. NEUTRON STAR - NEUTRON STAR COALESCENCES 

The most reliable source for the LIGO/VIRGO/GEO network is the coales- 
cences of neutron star/neutron star (NS/NS) binaries; this type of source 
has been studied theoretically in great detail over the last several years. 
The event rate is fairly well understood from observations of progenitor 
NS/NS systems in our own Galaxy: the distance to which one must look 
to see 3 events per year has been estimated to be ~ 200 Mpc [32, 33]. 
[Recent revisions to the pulsar distance scale indicate that this distance 
could be somewhat larger, perhaps ~ 400 Mpc [34]]. The range of initial 
LIGO interferometers for NS/NS binaries is ~ 25 Mpc, while that of the 
enhanced interferometers will be ~ 250 Mpc. There is a clean separation in 
the gravitational waveforms between an inspiral phase of the signal, which 
carries most of the detection signal-to-noise ratio and which is fairly well 
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understood theoretically, and a subsequent merger phase, which depends 
on the internal structure of the stars and which is poorly understood the- 
oretically [35]. Fortunately, understanding the details of the merger phase 
is not necessary for detecting the signals; for detection one will use the 
inspiral waves. 

Detection of gravitational waves from NS/NS binaries would test the 
fairly popular theory that NS/NS mergers are the source of gamma ray 
bursts, which are seen at a rate of roughly once per day by satellite gamma 
ray detectors, and whose origin has been a longstanding puzzle. The rate 
of burst detections of a few per day is roughly consistent with the inferred 
rate of NS/NS mergers out to cosmological distances, and amount of energy 
supplied by a NS/NS merger is (marginally) sufficient to produce the ob- 
served 7 ray energy fluxes for a source at a cosmological distance. In the last 
year new evidence has been uncovered that gamma ray bursts are cosmo- 
logical in origin rather than Galactic, which favors the NS/NS coalescence 
hypothesis: optical counterparts of some particular bursts have been seen 
at cosmological distances [36]. On the other hand, while the light curves 
and spectra of the bursts seem to be well fit by a model of an expanding 
relativistic fireball, there are theoretical difficulties in understanding how 
NS/NS mergers (or NS/BH or BH/BH mergers) could give rise to such fire- 
balls with sufficiently low baryon contamination. [Too many baryons in the 
fireball would convert too much of the fireball's energy into kinetic energy 
of the baryons instead of into radiation] . 

During the inspiral of the two stars, general relativistic effects are large 
as the orbital velocity is about ~ 20% of the speed of light at the point 
where most of the detection signal-to- noise ratio is obtained. The signal is 
very clean however: to a very good approximation, the two stars can be 
treated as point masses with intrinsic spins. The gravitational waveform is 
then determined by the relativistic gravitational interaction between the 
two point masses; all of the other complicated physical effects that occur 
(heating and tidal distortion of the stars, interactions involving the stars' 
magnetic fields etc.) have a negligible effect on the gravitational waveforms 
in the inspiral phase [37]. 

This cleanliness means that the characteristics of the waveforms can be 
predicted in advance, and that the technique of matched filtering [1] can be 
used to detect the waves. By cross-correlating the outputs of the various in- 
terferometers in the network with theoretical template waveforms, one can 
both detect the waves and also measure the binaries parameters: masses, 
spins, orbital elements, sky location, distance etc. Very high accuracy tem- 
plates are not required to detect the waves: currently available, post- 2- 
Newtonian templates will be adequate in most cases [38]. By contrast, in 
order to extract the parameters of the binary with acceptably small system- 
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atic errors, very high accuracy templates will be required; such templates 
are under construction [39] . It is possible that approximate post-Newtonian 
waveforms may be substantially improved using the method of Pade ap- 
proximants [40]. Much theoretical work has been done on how accurately 
the parameters of the binary can be measured; however, the parameter- 
extraction accuracies are currently only known to within factors of order 
2. 

There are several payoffs of measuring the inspiral waves from coalesc- 
ing binaries. One will be able to experimentally measure for the first time 
various non-linear aspects of general relativity such as gravitational wave 
tails [41, 42] and possibly in some cases the dragging of inertial frames [43]. 
It will be possible to constrain non general-relativistic theories of gravity, 
for example one will be able to improve on solar system bounds on the 
dimensionless parameter to of Brans-Dicke theory [44] , and one will obtain 
constraints on the "graviton mass" [45]. One could identify the sources of 
gamma ray bursts by comparing their arrival times to those of gravitational 
wave bursts. If the event rate is large enough, one will be able to accumu- 
late statistics on the distances, approximate sky locations, and masses of 
many (~ 100) coalescences. This information will directly probe the dis- 
tribution of Galaxies on scales 400 Mpc <^ D <, 4 Gpc with a resolution of 
~ 200 Mpc [41]; the distribution of Galaxies on such large scales has not 
yet been probed in redshift surveys. In addition, one will be able to at- 
tempt to measure the three cosmological parameters: the Hubble constant, 
the cosmological constant, and the closure parameter fi. Monte Carlo sim- 
ulations suggest that with one year's observations with advanced LIGO 
interferometers, one will be able to measure the Hubble constant to ~ 1% 
and the other two parameters with fractional accuracies of a few tens of 
percent [46]. [Of course, it is likely that these parameters will already have 
been accurately measured by the time advanced interferometer sensitivity 
levels are achieved, by the MAP and PLANCK satellite measurements of 
the CMBR [47]]. 

Turn, now, to the gravitational waves from the merger stage of coalesc- 
ing binaries. These will carry information about the internal structure of 
the neutron stars and potentially about the equation of state of nuclear 
matter at high densities. Efforts are underway to understand the dynamics 
of the merger process and the features of the emitted waves using numeri- 
cal simulations. So far, simulations of the merger have been limited to the 
Newtonian or post-l-Newtonian approximations [48]. There are indications 
that it might be necessary to simulate the merger numerically in full gen- 
eral relativity; efforts are underway in this direction. The merger waveforms 
will have most of their power at high frequencies of several hundred Hz or 
higher. Unfortunately, broadband laser interferometers have poor perfor- 
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Figure 4- A summary of the various phases of the waves (inspiral, merger, and ringdown) 
for the three different types of coalescing binary, and the usefulness of each phase for each 
type of source. 

mance at these high frequencies (see Fig. 1 above). However, it might be 
possible to measure the merger waves with specialized narrowband interfer- 
ometers which are designed to have enhanced sensitivity near some chosen 
frequency and worsened sensitivities elsewhere [49], or with the next gen- 
eration of resonant mass antennae (TIG As or spheres). 

6.2. NS/BH AND BH/BH BINARIES 

The other types of coalescing compact binary that might be measured are 
NS/BH and BH/BH binaries, as illustrated in Fig. 4. For the NS/BH case 
the signal can be divided into inspiral and merger phases as before; for the 
BH/BH case the signal can be divided into three phases, inspiral, merger 
and ringdown [50]. 

There are several fundamental differences between the NS/BH and BH/BH 
cases, and the NS/NS case: 

First, the event rates for NS/BH and BH/BH coalescences are far less 
certain than NS/NS coalescences, since so far we have no observed progeni- 
tor systems. Arguments based on progenitor evolution scenarios suggest an 
event rate rate very roughly ~ 3yr _1 within 200 Mpc, i.e., the same rate 
as for NS/NS binaries, in both the NB/BH and BH/BH cases [32]. [A re- 
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Figure 5. A plot of the wave strength y/Eh c g for waves from a binary black hole 
coalescence of total mass WOMq at redshift z = 1, together with the noise levels VEh lrnB 
for initial and advanced LIGO interferometers (adapted from Ref. [50]). The energy 
spectra for the merger waves and the ringdown waves are taken from the crude models of 
Ref. [50] (10% of total mass energy radiated in merger waves and 3% in ringdown waves). 
Notice that the inspiral waves are not visible with the initial LIGO interferometers, and 
barely visible with the advanced. For this type of source, it will be necessary to try to 
detect the waves using the merger or ringdown waves directly. 

cent estimate by Brown and Bethe reduces this distance to 100 Mpc for the 
NS/BH case [51]]. There are three classes of BH/BH binaries: field binaries, 
binaries formed by capture processes in galactic centers, and capture bina- 
ries in globular clusters. Siggurdson and Hernquist have argued that there 
is generically at least one BH/BH event per core collapse globular cluster 
[52]; this yields ~ 3 BH/BH events per year within a distance of 600 Mpc 
using the extrapolation method of Sec. 3.1 of Ref. [32]. For field BH/BH 
binaries, recent estimates of the coalescence rate by experts in binary evo- 
lution theory range from ~ 10 -8 yr _1 to ~ 10 -6 yr _1 in our Galaxy [53, 54], 
to completely negligible [55]. There are large uncertainties associated with 
these theoretical estimates of the coalescence rates. 

Second, for BH/BH coalescences, the merger waves will bring infor- 
mation about the dynamics of general relativity in a highly dynamical, 
highly nonlinear, highly non-spherical regime that has never before been 
probed experimentally. A major effort is underway in the numerical rela- 
tivity community to simulate the mergers of black holes, driven in part by 
the possibility of performing comparisons with gravitational wave data; see 
the article by Ed Seidel in this proceedings. The overall energy and the 
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amount of information carried by the merger waves are still very uncertain. 

Third, the inspiral waves are stronger for NS/BH and BH/BH sources 
than for NS/NS sources, simply because of their greater mass. The ini- 
tial LIGO interferometers can reach to ~ 40Mpc for a BH/BH binary of 
total mass M = 5M Q and to ~ 300 Mpc for M = 5OM [50]; enhanced 
interferometers can see roughly ten times further and advanced interferom- 
eters about 30 times further. Thus, BH/BH coalescences might be detected 
early by the LIGO/VIRGO network, and may be seen before NS/NS co- 
alescences. Note also that if BH/BH events are seen, the most commonly 
detected events will be the largest mass systems, since the increase in signal- 
to-noise ratio with increasing mass very likely more than compensates for 
the decrease in the mass function with increasing mass [56]. 

Fourth, the frequency at which the inpiral ends, i.e., at which the 
spinning-point-mass approximation fails, decreases as the total mass of the 
system is increased. For sufficiently massive systems, the regime of optimum 
sensitivity for LIGO (/ ~ 100 Hz) overlaps with the complicated merger 
waves and not with the simple inspiral waves; see Fig. 5. For such systems 
one needs to detect the events using the merger or ringdown waves; this 
implies an extra importance for numerical relativity simulations [50]. For 
black hole binaries of total mass ~ 2OM say (a conventional value) , events 
can be detected using inspiral waves; however a large part of the signal- 
to-noise ratio is accumulated in a very relativistic regime 6M r <^ 12M 
in which post-Newtonian detection templates are possibly inadequate. New 
numerical techniques are under development for calculating waveforms in 
this so-called "Intermediate Binary Black Hole" (IBBH) regime [57]. 

7. Periodic sources 

As is the case for supernovae, the expected event detection rate for periodic 
sources is highly uncertain due to our ignorance as to the wave's sources. 
The primary expected type of source are spinning neutron stars. Such stars 
can radiate in several ways [58]: 

First, if the axis of rotation is displaced from the corresponding principal 
axis of the moment of inertia tensor by a small "wobble angle" 8 W , then the 
star radiates at the precessional sideband of the rotation frequency with 
amplitude proportional to 9 w e p , where e p = {I xx — I Z z)/Izz is the poloidal 
eccentricity [59]. (We assume the axis of rotation is the z axis). However, 
the wobble angle decreases exponentially due to gravitational wave emission 
on the short timescale [60] 
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where / ro t is the rotational frequency. [See also Ref. [61] for another dissipa- 
tive mechanism which could also quickly damp away precessional motion.] 
Thus, this mechanism for producing waves is not likely to be detected, un- 
less possibly if the neutron star is accreting and if the accretion continuously 
drives 6 W away from zero, as has been suggested by Schutz [58]. 

Second, it is possible for the solid crust of the neutron star to support 
some non-axisymmetry. If the star is not axisymmetric, it will radiate at 
twice the rotation frequency with an amplitude proportional to the equato- 
rial ellipticity e e = (I xx — I yy )/I zz [59]. The effective strength of the waves 
is given by 

Here T G b s is the observation time and r is the distance to the source. For 
advanced LIGO interferometers it follows that the minimum detectable 
eccentricity is of order 

With specialized narrowband interferometers [49] or resonant mass detec- 
tors [12], this lower bound might be improved by two orders of magnitude 
[4]- " 

We are very ignorant as to the likely size of e e in spinning neutron stars. 
Theoretical estimates of the breaking strain of the star's crust yield the up- 
per limit e e ^ 10~ 5 [62]. The observed spin down rate of pulsars sets an 
upper limit on the eccentricity (obtained by assuming the entire spin down 
is due to gravitational wave emission); this limit is of order 10 -8 for mil- 
lisecond pulsars [1]. However, it does not follow that e e should be this small 
for other populations of neutron stars, whose typical evolutionary histories 
would be expected to be very different from those of millisecond pulsars. 
Several different plausible mechanisms for generating equatorial eccentric- 
ity have been suggested, for example the anisotropic pressure caused by 
magnetic fields [2, 63]. 

In the case of accreting systems (low mass X ray binaries), Bildsten has 
recently suggested that temperature gradients due to non-uniform accre- 
tion over the surface of the star lead (via temperature-dependent nuclear 
reactions) to density variations in the crust and thence to a non-zero e e ; the 
resulting estimated values of e e are of order 1CT 7 [64]. Bildsten also sug- 
gested an explanation for the observational puzzle that the rotation rate 
of several observed such systems cluster around 300 Hz: gravitational wave 
emission is preventing these sources from being spun up any further, i.e., all 
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the angular momentum being accreted is being radiated into gravitational 
waves. If this turns out to be true, several known systems such as SCO-X1 
would be detectable by advanced LIGO interferometers [64]. 

A third mechanism by which spinning neutron stars can emit gravi- 
tational waves is via the existence of unstable modes of vibration of the 
star. The /-mode of neutron stars will be unstable to a radiation-reaction 
driven instability called the Chandrasekhar-Friedman-Schutz instability at 
sufficiently rapid rotation rates (<; 0.9f2breakup)- ^ nas been shown that this 
instability is counteracted by the dissipative effects of viscosity, except when 
the star's temperature lies in the narrow band 10 9 K <> T <> 1Q W K [65]. 
Neutron stars are this hot only in the first few years after their formation, 
and newly born stars are typically not likely to be born spinning rapidly 
enough for the instability to operate. Thus, /-mode, radiation-reaction 
driven instabilities are not likely to be detected, except possibly in hot, 
accreting systems [66]. In addition, there are viscosity-driven "bar mode" 
instabilities which occur only for certain equations of state and for very 
rapidly rotating (P ~ 1ms) stars [63]; again, such instabilities are unlikely 
to be observable for the same reasons as for the CFS instability. 

Recently Anderson [67] , and also Friedman and Morsink [68] , discovered 
that neutron star r modes suffer from an analogous radiation-reaction- 
driven instability. Calculations indicate that the lowest r mode is unstable 
for temperatures T in the roughly the same range as above, when f2 <^ 20 Hz 
or so [69]. Thus, newly born neutron stars should be copious emitters of 
gravitational waves, and that they should always spin down from their ini- 
tial rotation rate to ~ 20 Hz on the timescale of a year after their formation 
[69]. Such strong sources should likely be visible to the VIRGO cluster of 
galaxies with initial LIGO interferometers, yielding an event rate of several 
per year. It is important to firm up the initial, crude estimates of wave 
strengths from these sources to verify this scenario. 

Finally, one should note that searching for neutron stars or any other 
type of periodic source is enormously computationally demanding due to 
the necessity of correcting for Doppler shifts due to the Earths motion 
and rotation, for many, many individual directions on the sky [70]. All 
sky searches will be computation limited: the minimum detectable signal 
strength with a Teraflop computer will be a factor of a few times larger than 
the ideal minimum detectable signal strength one would obtain with infi- 
nite computing power. Efforts are underway to devise sophisticated search 
algorithms to minimize the loss in detectable signal strength [70, 58]. 
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8. Low Frequency Sources: Supermassive Black Holes 

Consider now the low frequency band 1CP 4 Hz ^ / ^ 0.1 Hz, the domain of 
the proposed spacebased interferometer LISA. There is a guaranteed class 
of sources in this frequency band: binary stars, including common envelope 
binaries, close white dwarf binaries, contact binaries, cataclysmic variables, 
etc. For many binary, one should be able to combine optical and gravita- 
tional wave information to solve for all the orbital elements of each binary 
and to test general relativity. At low frequencies, there will be so many white 
dwarf binaries that they will be experimentally indistinguishable and will 
constitute a source of background noise. Also, the shortest period NS/NS 
binary in the Galaxy should have a period of about ~ 10 5 yrs, from the 
estimates of NS/NS coalescence rates. Such a binary would be visible to 
LISA with a signal-to-noise ratio of ~ 500 [2]. 

LISA should also be able to see the formations and/or coalescences of su- 
permassive black holes (SMBHs) throughout the observable Universe with 
signal to noise ratios of hundreds or thousands. There is some observational 
evidence for SMBH binaries: wiggles in the radio jet of QSO 1928+738 have 
been attributed to the orbital motion of a SMBH binary [71], as have time 
variations in quasar luminosities [72] and in emission line redshifts [73]. 
The overall event rate is uncertain, but could be large 1/yr), especially 
if the hierarchical scenario for structure formation is correct [74] . Detecting 
coalescences of supermassive black holes would allow high precision tests 
of general relativity. 

From the waves from a small black hole or compact object spiraling into 
a larger black hole, one can in principle obtain detailed information about 
the spacetime geometry of the larger black hole (its multipole moments 
[75] ) , and test the no-hair theorem of general relativity [4] . As a foundation 
for extracting this information, theorists must understand how to calculate 
the motion of a test particle in a Kerr geometry under the influence of 
gravitational radiation reaction. Some progress has been made recently in 
this direction: a general expression for the radiation reaction force on a test 
particle in any background spacetime has been derived by three indepen- 
dent methods, all of which give the same answer [76]. Ori [77] has proposed, 
without proof, a practical calculational scheme to calculate the evolution 
of orbits: it would be useful to verify that Ori's scheme can be derived from 
the general Mino-Quinn-Wald equation of motion [76] . 

9. Sources of Stochastic Gravitational Waves 

I will not attempt to give a review of stochastic sources here, but instead 
refer the reader to the recent detailed and comprehensive review article by 
Allen [5]. 
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10. Conclusion 

With the kilometer scale interferometer network about to come online in 
the next few years, this is an exciting time for gravitational wave astron- 
omy. Almost certainly we have not anticipated all the detectable sources of 
gravitational waves that Nature produces in the real Universe. Hopefully 
she will bring us some surprises. 
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